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Abstract 

We  study  the  response  to  perturbation  of  non-Poisson  dichotomous  fluctuations  that 
generate  super-diffusion.  We  adopt  the  Liouville  perspective  and  with  it  a  quantum-like 
approach  based  on  splitting  the  density  distribution  into  a  symmetric  and  an  anti-symmetric 
component.  To  accomodate  the  equilibrium  condition  behind  the  stationary  correlation 
function,  we  study  the  time  evolution  of  the  anti-symmetric  component,  while  keeping  the 
symmetric  component  at  equilibrium.  For  any  realistic  form  of  the  perturbed  distribution 
density  we  expect  a  breakdown  of  the  Onsager  principle,  namely,  of  the  property  that  the 
subsequent  regression  of  the  perturbation  to  equilibrium  is  identical  to  the  correspond¬ 
ing  equilibrium  correlation  function.  We  find  the  directions  to  follow  for  the  calculation  of 
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higher-order  correlation  functions,  an  unsettled  problem,  which  has  been  addressed  in  the  past 
by  means  of  approximations  yielding  quite  different  physical  effects. 

©  2004  Elsevier  B.V.  All  rights  reserved. 
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Keywords:  Stochastic  processes;  Non-Poisson  processes;  Liouville  and  Liouville-like  equations; 
Correlation  function;  Regression  to  equilibrium 


1.  Introduction 

Complex  physical  systems  typically  have  both  nonlinear  dynamical  and 
stochastic  components,  with  neither  one  dominating.  The  response  of  such 
systems  to  external  perturbations  determines  the  measurable  characteristics  of 
phenomena  from  the  beating  of  the  human  heart  to  the  relaxation  of  stressed 
polymers.  What  distinguishes  such  complex  phenomena  from  processes 
successfully  studied  using  equilibrium  statistical  mechanics  is  how  these  systems 
internalize  and  respond  to  environmental  changes.  Consequently,  it  is  of 
broad  interest  to  determine  which  of  the  prescriptions  from  equilibrium 
statistical  physics  is  still  applicable  to  complex  dynamical  phenomena  and 
which  are  not.  Herein  we  address  the  breakdown  of  one  of  these  fundamental 
relations,  that  is,  the  Onsager  Principle.  The  authors  of  Ref.  [1]  discussed  how 
to  make  the  generalized  master  equation  (GME)  compatible  with  the  Onsager 
principle,  leaving  open,  however,  the  practical  problem  of  the  response  to 
realistic  perturbations  that,  as  we  shall  see,  cause  the  breakdown  of  this  important 
property. 

It  is  worthwhile  to  briefly  review  this  issue  in  the  standard  cases.  In  the  case  of 
ordinary  statistical  mechanics  an  exhaustive  treatment  of  the  relaxation  of 
perturbations  to  equilibrium  can  be  found  in  Ref.  [2].  Let  us  consider  as  a  prototype 
of  ordinary  statistical  mechanics  the  case  when  the  stochastic  variable  under  study 
( l(t)  is  described  by  the  linear  Langevin  equation 


—  =  -y£0)  +  r\(t)  , 


(1) 


where  the  random  driving  force  rj(t)  is  white  noise.  Let  us  imagine  that  £(Y)  is  the 
velocity  of  a  particle  with  unit  mass  and  a  given  electrical  charge.  Furthermore,  we 
assume  that  this  system  reaches  the  condition  of  equilibrium,  and  at  a  given  time 
t  =  0,  we  apply  an  electrical  field  E(t).  The  external  field  E(t)  is  an  arbitrary  function 
of  time,  fitting  the  condition  that  E(t)  —  0,  for  t<  0.  The  adoption  of  linear  response 
theory  yields  the  prescription  for  the  mean  response  of  the  system  to  the  external 


field 


(2) 
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where  is  the  equilibrium  correlation  function  of  t.  In  this  familiar  case,  of 
dissipative  Brownian  motion,  the  correlation  function  of  the  particle  velocity  is  the 
exponential  exp(— yt). 

There  are  two  limiting  cases  of  time  dependence  of  the  perturbation:  (a)  the 
electric  field  is  proportional  to  the  Heaviside  step  function,  E(t )  =  K0(t),  and  is 
therefore  a  constant  field  after  it  is  turned  on  at  t  =  0;  (b)  the  electric  field  is 
proportional  to  the  Dirac  delta  function,  E(t)  =  K5(t),  and  is  consequently  an  initial 
pulse  that  perturbs  the  particle  velocity.  In  these  two  limiting  cases  we  obtain  for  the 
velocity  of  the  Brownian  particle 

m)  =  K  f  ^dt'  (3) 

Jo 

and 


m)  =  K^(t) ,  (4) 

respectively.  These  two  limiting  cases  show  that  in  the  case  of  ordinary  statistical 
mechanics  the  system’s  response  to  an  external  perturbation  is  expressed  in  terms  of 
the  unperturbed  correlation  function.  We  shall  refer  to  the  conditions  of  Eqs.  (3)  and 
(4)  as  the  Green-Kubo  relation  and  the  Onsager  relation,  respectively. 

The  search  for  a  dynamical  derivation  of  anomalous  diffusion,  that  is,  where  the 
mean  square  value  of  the  dynamic  variable  is  not  linear  in  time,  has  been  a  subject  of 
great  interest  in  recent  years.  There  are  two  main  theoretical  perspectives  on  how  to 
explain  the  origin  of  anomalous  diffusion.  The  first  perspective  is  based  on  the 
assumption  that  there  are  unpredictable  events,  that  the  occurrence  of  these  events 
obey  non-Poisson  statistics,  and  is  related  to  the  pioneering  paper  by  Montroll  and 
Weiss  [3].  For  more  recent  research  work  on  this  subject  we  refer  the  reader  to  the 
excellent  review  papers  of  Refs.  [4,5].  The  other  perspective  rests  on  the  assumption 
that  single  diffusion  trajectories  have  an  infinite  memory.  The  prototype  of  the  latter 
perspective  is  the  concept  of  fractional  Brownian  motion  introduced  by  Mandelbrot 
[6],  A  problem  worthy  of  investigation  is  whether  or  not,  in  the  case  of  anomalous 
diffusion,  the  response  to  external  perturbation  departs  from  the  predictions  of  Eqs. 
(3)  and  (4).  In  the  last  few  years,  this  problem  has  been  addressed  by  some 
investigators  [7-12].  These  authors  have  discussed  the  Green-Kubo  relation  of  Eq. 
(3).  Notice  that  in  the  special  case  of  ordinary  statistical  mechanics  this  relation  can 
also  be  written  in  the  following  form 

M0>  =  t-^27  {x2(t))0  .  (5) 

To  understand  how  to  derive  this  equation,  originally  proposed  by  Bouchaud  and 
George  [13],  we  refer  to  the  following  equation  of  motion: 


dx 

d? 


=  m . 


(6) 


Since  £( t )  is  a  fluctuating  velocity,  it  generates  spatial  diffusion  and  we  denote  by  x(t) 
the  position  of  the  corresponding  diffusing  particle.  The  external  field  affects  the 
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velocity  fluctuation  and,  consequently,  the  diffusion  process  generated  by  these 
fluctuations.  In  the  absence  of  perturbations,  the  second  moment  of  the  diffusing 
particle,  {x(t)2)0,  obeys  the  prescription 


(7) 


It  is  straightforward  to  prove  that  Eq.  (3)  yields  Eq.  (5).  This  is  done  by  considering 
the  mean  value  of  Eq.  (6) 


(8) 


The  time  integration  of  the  left-hand  term  of  this  equation  yields  the  left-hand  term 
of  Eq.  (5)  and  the  time  integration  of  the  right-hand  term  of  Eq.  (8),  using  Eqs.  (3) 
and  (7),  yields  the  right-hand  term  of  Eq.  (5). 

The  relation  of  Eq.  (5)  is  denoted  as  generalized  Einstein  relation ,  because  it  might 
hold  true  also  when  the  equilibrium  correlation  function  does  not  exist  [10]. 
However,  in  the  case  of  ordinary  statistical  mechanics  Eq.  (5)  becomes  equivalent  to 
the  Green-Kubo  property.  In  this  generalized  sense  we  can  state  that  the  authors  of 
Refs.  [7-12]  studied  the  Green-Kubo  relation  of  Eq.  (3),  all  of  them  but  the  authors 
of  Ref.  [7],  devoting  their  attention  to  the  sub-diffusional  case,  that  is,  the  mean- 
square  average  displacement  is  proportional  to  f  with  y<  1.  It  is  worth  mentioning 
that  there  exists  experimental  evidence  in  favor  of  this  relation  [14,15]. 

Herein  we  focus  our  attention  on  the  Onsager  relation  of  Eq.  (4).  The  only 
earlier  work  on  how  to  extend  the  Onsager  principle  to  the  case  of  non-Poisson 
statistics,  known  to  us,  is  that  of  Ref.  [1].  In  this  paper  the  problem  is  addressed 
analyzing  a  virtually  infinite  symbolic  sequence,  of  +’s  and  — ’s,  corresponding  to  the 
two  values  of  £(?),  which,  in  fact,  is  supposed  to  be  a  dichotomous  variable.  They 
made  the  ergodic  and  stationary  assumption  and  proved  that  this  assumption  yields 
the  Onsager’s  regression  principle.  Adopting  the  same  view  with  a  Gibbs  set  of 
infinitely  many  independent  sequences,  the  prescription  of  Ref.  [1]  would  work  as 
follows.  We  set  the  preparation  at  t  =  —  oo.  This  ensures  the  equilibrium  condition, 
if,  as  assumed  in  Ref.  [1],  the  non-Poisson  statistics  is  compatible  with  ergodicity. 
Then  at  a  given  time  t>0,  we  select  only  the  systems  where  ^  —  +.  From  this  time 
on,  we  follow  the  time  evolution  of  these  systems  and  we  prove  that  (£(?)}  is 
proportional  to  the  correlation  function  </>;(;),  implying  the  validity  of  Onsager’s 
principle. 

This  conclusion  is  correct,  but  the  choice  done  to  prepare  the  initial  out-of¬ 
equilibrium  condition  state  implies  that  we  have  a  detailed  information  on  the  single 
constituents  of  the  Gibbs  ensemble,  an  ideal  condition  difficult,  if  not  impossible,  to 
realize  in  practice.  Here  we  plan  to  address  the  same  problem  from  a  perspective 
compatible,  in  principle,  with  experiment,  by  adopting  a  Liouville-like  approach. 
This  approach  will  allow  us  to  establish  some  general  conclusions  about  the  issue 
raised  by  Bologna  et  al.  [16]  about  a  possible  conflict  between  the  density  and 
trajectory  picture.  For  this  reason,  we  shall  illustrate  the  rules  for  the  calculation  of 
the  correlation  function  $<_=(?),  with  a  prescription  that  can  be  in  principle  extended  to 
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correlation  functions  of  any  order.  We  expect  that  these  prescriptions  might  lead  to  a 
successful  evaluation  of  the  fourth-order  correlation  function,  which  has  been 
studied  so  far  by  means  of  a  factorization  assumption  which  is  violated  by  the  non- 
Poisson  statistics. 


2.  An  idealized  model  of  intermittent  randomness  and  the  corresponding  density 
equation 

As  done  in  Ref.  [16],  let  us  focus  on  the  following  dynamical  system.  Let  us 
consider  a  variable  y  moving  within  the  interval  /  =  [0, 2].  The  interval  is 
defined  over  an  overdamped  potential  V,  with  a  cusp-like  minimum  located  at 
y  =  1.  If  the  initial  condition  of  the  particle  is  y( 0)  >  1 ,  the  particle  moves  from  the 
right  to  the  left  toward  the  potential  minimum.  If  the  initial  condition  is  y(0)<l, 
then  the  motion  of  the  particle  toward  the  potential  minimum  takes  place 
from  the  left  to  the  right.  When  the  particle  reaches  the  potential  bottom  it  is 
injected  to  an  initial  condition,  different  from  y  =  1,  chosen  in  a  random 
manner.  We  thus  realize  a  mixture  of  randomness  and  slow  deterministic 
dynamics.  The  left  and  right  portions  of  the  potential  V (y)  correspond  to  the 
laminar  regions  of  turbulent  dynamics,  while  randomness  is  concentrated  at  y  =  1. 
In  other  words,  this  is  an  idealization  of  the  map  used  by  Zumofen  and  Klafter  [17], 
which  does  not  affect  the  long-time  dynamics  of  the  process,  yielding  only  the  benefit 
of  a  clear  distinction  between  random  and  deterministic  dynamics.  Note  that  the 
waiting  time  distribution  in  the  two  laminar  phases  of  the  reduced  form  has  the  same 
time  asymptotic  form  as 


(9) 


We  select  this  form  as  the  simplest  possible  way  to  ensure  the  normalization 
condition 


(10) 


We  note  that  Eq.  (10)  implies  )i>  1.  The  condition  /r>2  corresponds  to  the  existence 
of  a  finite  mean  sojourn  time,  and,  thus,  to  the  possibility  itself  of  defining  the 
stationary  correlation  function  of  the  fluctuation  c,  which,  with  the  choice  of  Eq.  (9) 
reads  [18] 


From  within  the  perspective  of  a  single  trajectory  this  dynamical  model  has  the 
form 


y  =  'W  i  -  >’)/  -  &(y  -  m  -  yf]  +  - —  S(y-  l) . 


^ random 


(12) 
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The  function  0(x)  is  the  ordinary  Heaviside  step  function,  Av(t)  is  a  random  function 
of  time  that  can  achieve  any  value  on  the  interval  [-1,4-1],  and  r random  is  the  injection 
time  that  must  fulfill  the  condition  of  being  infinitely  smaller  than  the  time  of  sojourn 
in  one  of  the  laminar  phases.  Note  that  z  is  a  real  number  fitting  the  condition  z>  1. 
In  fact,  the  equality 


(h  -  1) 


(13) 


relates  the  dynamics  of  Eq.  (12)  to  the  distribution  Eq.  (9).  The  Poisson  condition  is 
recovered  in  the  limit  z  — >  1,  namely,  in  the  limit  /(  — »■  oo.  Thus,  in  a  sense,  the  whole 
region  z>l  (p<oo)  corresponds  to  anomalous  statistical  mechanics.  However,  the 
deviation  from  normal  statistical  mechanics  is  especially  evident  when  z>1.5,  a 
condition  implying  that  the  second  moment  diverges.  In  the  case  z  >  2  the  departure 
from  ordinary  statistical  mechanics  becomes  even  more  dramatic,  due  to  the  fact  that 
the  first  moment  also  diverges  and,  as  we  shall  see  in  this  Section,  the  process 
becomes  non-ergodic. 

Let  us  now  move  to  the  density  picture,  namely,  to  a  formulation  of  Eq.  (12)  from 
within  the  Gibbs  perspective.  The  form  of  this  equation  of  evolution  for  the 
probability  density  is 

^ p(y,  t )  =  -A  ^  {[0(1  -  >’)/  -  &(y  -  i)(2  -  yf]p(y,  0}  +  C(t) ,  (14) 

where 


C(t)  =  kp(\,t)  . 


(15) 


It  is  important  to  stress  that  we  are  forced  to  set  the  equality  of  Eq.  (15)  to  fulfill  the 
following  physical  conditions: 


d 

d  7 


//=[  0,2] 


p(y,  t)  d y 


ll=[  0,2] 


0 

07 


p(v,  7)  dy  =  0  , 


(16) 


which,  in  turn,  ensures  the  conservation  of  probability.  We  assume  the  ordinary 
normalization  condition 


/  p(y,  t)  dv  =  1  ,  (17) 

.//— [0,2] 

which  is  kept  constant  in  time,  as  a  consequence  of  Eq.  (16).  It  is  evident  that  the 
inhomogeneous  term  C(7)  corresponds  to  the  action  of  the  stochastic  term,  namely, 
the  second  term  on  the  right-hand  side  of  Eq.  (12). 

It  is  important  to  point  out  that  our  dynamic  perspective  allowed  us  to  describe 
the  intermittent  process  through  the  Liouville-like  equation 


9pQ’,  t) 

07 


=  @p(y,  t ) , 


(18) 


where  y  denotes  a  continuous  variable  moving  either  in  the  right  or  in  the  left 
laminar  region,  with  £  getting  the  values  IE  or  —IE,  correspondingly,  and  the 
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operator  M  reading 


=  -2  A  [0(1  -  y)f  -  &(y  -  1 )( 2  -  yf]  +  k  J*  dy«5(y  -  1)  . 


(19) 


This  operator  departs  from  the  conventional  form  of  a  differential  operator, 
since  the  last  term  corresponds  to  the  unusual  role  of  an  injection  process, 
which  is  random  rather  than  being  deterministic.  In  the  ordinary  Fokker-Planck 
approach  the  role  of  the  stochastic  force  is  played  by  a  second-order  derivative, 
which  is  not  as  unusual  as  the  integral  operator  of  Eq.  (19).  Here  the 
role  of  randomness  is  played  by  the  back-injection  process,  which,  from 
within  the  density  perspective  is  described  by  an  operator  that  selects 
from  all  possible  values  p(y,  t),  the  specific  value  of  p(y,  t)  at  y  =  1 .  The  idealization 
that  we  have  adopted,  of  reducing  the  size  of  the  chaotic  region  to  zero, 
with  the  choice  of  the  process  of  back  injection  located  at  y  =  1,  makes  it 
possible  for  us  to  use  the  continuous  time  representation  and  the  equation  of  motion 
Eq.  (18)  rather  than  the  conventional  Frobenius-Perron  representation.  This 
representation  will  allow  us  to  obtain  analytical  results.  However,  the  same  physical 
conclusions  would  be  reached,  albeit  with  more  extensive  algebra,  using  the 
conventional  maps  and  the  Frobenius-Perron  procedure  described  in  the  recent 
book  by  Driebe  [19]. 

We  note  that  the  equilibrium  probability  density  solving  Eq.  (14)  is  given  by 


Po(y)  = 


2  —  z 
2 


0(1  -y)  0(y-l)' 

.  f-1  (2  -yf-\  ' 


(20) 


This  equilibrium  density  becomes  negative  for  z>2.  This  fact,  by  itself,  is  not  a  big 
problem.  As  we  shall  see  in  Section  5,  we  might  define  the  distribution  density  as 
—p0(y),  thereby  settling  this  problem.  There  is,  however,  a  deeper  reason  why  we 
cannot  rest  in  Eq.  (20),  when  z>  2.  This  is  the  fact  that  the  distribution  cannot  be 
normalized,  thereby  signaling  the  important  fact  that  for  r  >  2  there  no  longer  exists 
an  invariant  distribution.  The  lack  of  an  invariant  distribution  accounts  for  the  non- 
ergodicity  in  the  fluorescence  of  single  nanocrystals,  modeled  as  a  dichotomous 
process,  recently  pointed  out  by  Brokmann  et  al.  [20]. 

For  the  purposes  of  calculation  in  the  next  few  sections,  it  is  convenient  to 
split  the  density  p(y,  t)  into  a  symmetric  and  an  anti- symmetric  part  with 
respect  to  y  =  1 , 


pi}’,  *)  =  Psiy >  0  +  Pa  (y,  0  ■  (21) 

This  separation  based  on  symmetry  yields  the  following  two  equations  from  Eq.  (14): 
~tPs(y,  t)  =  -k&(  1  -  y)  \fPsiy,  0]  +  k&(y  -  1)  ^  [(2  -  y)zps(y,  r)]  +  C(f) 


(22) 
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and 

^PaLv,  t)=  -W(i-  y)  ^  {ypA(y,  0] 

+  W(y  -  1)  ^  [(2  -  yfpA(y,  t)\  .  (23) 

We  note  that  the  anti-symmetric  part  of  the  density  is  driven  by  a  conventional 
differential  operator,  which  we  denote  by  f .  Thus,  we  rewrite  Eq.  (23)  as  follows: 

^PaLvJ)  =  rpA(y,t)  ,  (24) 

where 


r  =  -A&(  l 


■y)—  y+W(y 
dy 


l)-(2  -yf 
dy 


Then,  it  is  convenient  to  define 


Ps,a(}’=  Ut)  = 


1 

2 


(25) 


(26) 


which  implies  pA(y  =  1,  t)  =  0. 

The  reader  might  be  disturbed  by  the  fact  that  these  results  have  been  derived  by 
dealing  with  the  step  function  as  a  constant,  while  it  apparently  yields  a  significant 
discontinuity  at  y  =  1 . 

Actually,  as  we  shall  see  in  Section  5,  the  discontinuity  affects  only  the  anti¬ 
symmetric  part,  and  this  discontinuity  is  justified  by  the  fact  that  the  right  and  left 
portions  of  the  anti-symmetric  part  correspond  to  two  independent  dynamic 
processes  of  particles  and  holes. 

The  operator  with  the  unusual  form,  containing  C(t),  is  only  responsible  for  the 
time  evolution  of  the  symmetric  part  of  the  probability  density.  We  notice,  on  the 
other  hand,  that  any  physical  effect  producing  a  departure  of  C(t)  from  its 
equilibrium  value,  if  this  exists,  namely,  if  r<2,  implies  a  departure  from 
equilibrium.  A  stationary  correlation  function  can  be  evaluated,  as  we  shall  see  in 
the  next  few  sections,  using  only  Eq.  (23),  without  forcing  Eq.  (22)  to  depart  from  the 
equilibrium  condition. 

As  we  shall  see  in  Section  6,  the  evaluation  of  correlation  functions  of  order  higher 
than  the  second  cannot  be  done  without  producing  a  departure  of  C(t)  from  its 
equilibrium  value.  This  might  generate  the  impression  that  the  correlation  functions 
of  order  higher  than  the  second  cannot  be  evaluated  without  internal  inconsistencies, 
if  we  use  only  the  density  picture.  The  evaluation  of  these  higher-order  correlation 
functions  was  done  in  Ref.  [21],  by  using  a  procedure  based  on  the  time  evolution  of 
single  trajectories.  Actually,  as  we  shall  see  in  Section  6,  the  density  approach  should 
yield  the  same  result.  However,  we  think  that  deriving  this  result  using  only  the 
Liouville-like  equation  of  this  section  is  a  hard  task,  which  was  bypassed  in  the 
past  by  means  of  the  factorization  approximation  [16,25],  violated  by  the  non- 
Poisson  case. 
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3.  The  correlation  function  of  the  dichotomous  fluctuation  from  the  trajectory  picture 


Let  us  focus  our  attention  on  Eq.  (12)  and  consider  the  initial  condition  y0  e  [0, 1]. 
Then,  it  is  straightforward  to  prove  that  the  solution,  for  y<  1  is, 

y(t)  =  y0(  1  -  A(z  -  1  )f0-lt)~mz-l)  .  (27) 

From  (27)  and  imposing  the  condition  y(T)  =  1,  we  can  find  the  time  at  which  the 
trajectory  reaches  the  point  y  =  1,  which  is 


T  =  T(y0 )  = 


Kz  —  i)jo_1 


(28) 


Since  we  have  to  find  £,(t)  and  c(t)  —  ^iy(t)),  from  the  general  form  of  Eq.  (28)  we 
obtain: 


^  =  [0(1  -  yo)0(T(yo )  -  t)  -  0(yo  -  1)0(T(2  -  v0)  -  0] 

+OQ 

-  X  si8n 
;=0 

(29) 


X>* 

\k= 0  ) 


0lV'Tk-t)-0  Y,Tk-t 


Z^ 

VA-0 


Z^ 

VA-0 


where  the  time  increments  are  given  by 


rr,  ,  1  -  Jo  1  nn  *  1  -  (2  -  Jo) 

to  =  Hjo)  =  77: — r^rr  0(1  -  j0)  + 


1 


2(z  -  1)Jq 


2(r-  l)(2-v0> 


Z —  1 


0(Jo  -  1) 


i  -  [i  +  A  ,  ^  ,  i  -  [i  -  -Ur,)]  1 

A(z  -  1)[1  +  Ay (Ur1  j  A(Z  -  1  )[1  -  I .(r, )]  1 


0(  MU)  . 


(30) 


Then,  for  the  correlation  function  we  obtain  the  following  expression: 

mm 


wz 


([0(1  -  Jo)0(njo)  -  0  +  0(Jo  -  1)0(212  -  v0)  -  0]> 


+oo 


X(  si§n(Jo  -  1  )siSn 


;'=0 


(  X  Xk 


\  k=0 


\k=0 


\k=0 


(31) 


As  pointed  out  in  Section  2,  the  calculation  of  the  correlation  function  rests  on 
averaging  on  the  invariant  distribution  given  by  Eq.  (20).  As  a  consequence  of  this 
averaging,  the  second  term  in  Eq.  (31)  vanishes.  In  fact,  the  quantity  to  average  is 
anti-symmetric,  whereas  the  statistical  weight  is  symmetric. 
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It  is  possible  to  write  the  surviving  term  in  the  correlation  as 


fl+t(z- 


=  (2  -  2) 


=  (1  +  A(z  -  l)t)-v- 
=  (i + , 


-(2-z)/(z-l) 


(32) 


with 


(33) 


Since  we  focus  our  attention  on  0</?<l,  we  have  to  consider  ^<r<2.  Note  that 
the  region  \<z<\  does  not  produce  evident  signs  of  deviation  from  ordinary 
statistics.  However,  as  we  shall  see  in  Section  6,  the  higher-order  correlation 
functions  can  be  easily  evaluated  only  in  the  case  z—  1 ,  when  the  correlation 
function  becomes  identical  to  the  exponential  function  exp(— It),  where  the 
factorization  assumption  holds  true.  If  z>  1,  the  calculation  of  the  higher-order 
correlation  functions  from  within  the  Liouville-like  approach  becomes  a  very 
difficult  task.  Note  also  that  Eq.  (32)  becomes  identical  to  Eq.  (11)  after  setting  the 
condition 


2(7  -  1)  =  y  ■ 


(34) 


4.  The  correlation  function  of  the  dichotomous  fluctuation  from  the  density  picture 

The  result  of  the  preceding  section  is  reassuring,  since  it  establishes  that  the 
intermittent  model  we  are  using  generates  the  wanted  inverse  power  law  form  for  the 
correlation  function  of  the  dichotomous  variable  £(f).  In  the  present  section  we  show 
that  exactly  the  same  result  can  be  derived  from  the  adoption  of  the 
Frobenius-Perron  form  of  Eq.  (14). 

To  fix  ideas,  let  us  concisely  review  the  results  obtained  in  earlier  work  [22,23], 
which  refers  to  a  dynamic  system  that  can  be  identified  with  the  left  part  of  the  model 
of  Section  2.  A  particle  in  the  interval  [0, 1]  moves  toward  y  =  1  following  the 
prescription  y  =  ).yz  and  when  it  reaches  y  =  1  it  is  injected  backwards  at  a  random 
position  in  the  interval.  The  evolution  equation  obeyed  by  the  densities  defined  on 
this  interval  is  the  same  as  Eq.  (14),  with  C(t)  =  /-/>(  1 ,  t).  This  dynamic  problem  was 
addressed  in  Refs.  [22,23],  and  solved  using  the  method  of  characteristics  as  detailed 
in  Ref.  [24].  It  is  important  to  stress  that  this  approach  is  the  requisite  price  for 
adopting  the  idealized  version  of  intermittency.  The  adoption  of  the  more 
conventional  reduced  map  of  Ref.  [17]  would  have  made  it  possible  for  us  to  adopt 
the  elegant  prescriptions  of  Driebe  [19],  as  done  in  Ref.  [22]. 
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Let  us  remind  the  reader  that  the  solution  afforded  by  the  method  of 
characteristics,  in  the  case  of  this  simple  non-linear  equation  with  stochastic 
boundary  conditions  is 


piy,  t)  = 


MU) 


I o  (iz(U  -  0yz~') 


dC+p 


y 


,0 


1 


9z(yz  xi) 


(35) 


where 


</z(x)  =  (l  +A(z-  \)x)z/(z-])  •  (36) 

As  shown  later  in  this  section,  to  find  the  correlation  function  (£(Y)<i;(0)},  and  the 
corresponding  functions  of  higher  order  as  well,  using  only  densities,  we  have  to 
solve  Eqs.  (22)  and  (23),  which  are  equations  of  the  same  form  as  that  yielding  Eq. 
(35).  For  this  reason  we  adopt  the  method  of  characteristics  again.  To  do  the 
calculation  in  this  case,  it  is  convenient  to  adopt  a  frame  symmetric  with  respect  to 
y  =  1.  Then,  let  us  define  Y  =  y  —  1.  Using  the  new  variable  and  Eq.  (35),  we  find  for 
Eqs.  (22)  and  (23)  the  following  solution: 


Ps(Y’t) 


f"  4ps(°»t)  dT 

Jo  cJz((t  -  r)(l  -  \Y\rl) 


Us 


1  - 


i-m 


M((i-iTir'o]1/z 


,o 


i 


g£(\-\Y\Tli) 


(37) 


and 


pA(Y,  t)  =  pA  ^sign(T) 

1 

x - — i—  •  (3S) 

g((l-\Y\r't) 

Then,  the  solution  consists  of  two  terms:  (1)  the  former  is  an  even  term  and  is 
responsible  for  the  long-time  limit  of  the  distribution  evolution  and  (2)  the  latter  is 
an  odd  term  which  disappears  in  the  long-time  limit.  We  note  that  (2)  is  a  desirable 
property  because  Eq.  (23)  does  not  contain  the  injection  term  C(t)  and  the 
equilibrium  density  (20)  is  an  even  function,  independently  of  the  symmetry  of  the 
initial  distribution. 

Let  us  illustrate  the  rules  that  we  adopt  in  this  paper  to  evaluate  the  equilibrium 
correlation  function  (£(f«)£(bi-i)  ■  ■  ■  £(h)),  with  ^  •  •  •  ~^t\  ^0.  We  assume 

the  system  to  be  prepared  in  the  equilibrium  state p0(Y),  which  is  identified  with  the 
Let  state  |0).  We  apply  to  it  the  time  evolution  operator  exp (fflt),  with  the  operator  3k 
defined  by  Eq.  (18).  This  does  not  have  any  effect,  since  formally  |0)  is  an  eigenstate 
of  9k-,  with  vanishing  eigenvalue  a  =  0.  Thus,  i2?|0)  =  a|0)  =  0.  At  time  t  —  t\  we 
apply  to  |0)  the  sign  operator,  which  yields  1  if  T<0  and  —1  if  T>0.  The 
eigenvalues  of  this  operator  coincide  with  the  values  of  the  stochastic  fluctuation  £(t). 
For  the  sake  of  simplicity  we  assign  to  this  operator  the  same  symbol  £,  as  the 
stochastic  fluctuation. 


1 


I  Y\ 


[gz((l-\Y\y-lt)\ 


,0 
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The  application  of  £  to  |0>  generates  an  anti-symmetric  distribution  pA(Y).  Then 
we  let  the  system  evolve  from  time  t—t\  to  ti,  with  the  time  evolution  obeying  Eq. 
(23),  which  is  formally  given  by  exp(T(?2  —  h  ))pA(Y).  At  time  t  =  ti  we  apply  to  this 
distribution  the  operator  £  again,  with  the  effect  of  changing  it  in  the  symmetric 
distribution  ps(  Y).  At  this  stage  we  let  this  distribution  evolve  in  time,  according  to 
Eq.  (14),  namely,  as  exp(^(?3  —  t2))ps(Y),  and  so  on.  The  calculation  is  concluded 
with  an  average  over  the  equilibrium  distribution  p0(  Y). 

In  the  ordinary  quantum  mechanical  case  where  the  dynamic  operators  are  super¬ 
operator  associated  to  Hermitian  Hamiltonian,  these  rules  yield  results  equivalent  to 
the  ordinary  prescription,  with  the  observable  motion  expressed  by  means  of  the 
Heisenberg  representation. 

As  pointed  out  in  Section  2,  according  to  these  rules  the  two-time  correlation 
function  is  determined  by  the  anti-symmetric  part  of  the  probability  density  alone, 
and  its  explicit  expression  becomes: 


mm)  =  (2  -  z) 


i 


0  (1-T) 


1 


T=i-(i-y)/[gz((i-y)I-1o]1/z  6b((l  —  0 


d  v 


(39) 


The  integral  (39)  is  exactly  solvable  and  leads  to  the  expression 

mm)  =  a  +  A(Z  - 1  )?)-(2-d/b-i) ,  (40) 

which  is  the  same  result  as  that  found  in  Section  3,  using  trajectories  rather  than  the 
probability  densities.  In  a  similar  way,  it  is  possible  to  calculate  the  correlation 
{Y(t)Y( 0))  and  determine  that  its  temporal  behavior  is  an  inverse  power  law  with  the 
same  exponent  as  that  in  Eq.  (40). 


5.  Onsager  regression  to  equilibrium 


In  conclusion,  in  the  two  preceding  sections  we  have  established  that  the  Liouville- 
like  representation  of  Eq.  (18)  yields,  as  expected,  the  equilibrium  correlation 
function  of  Eq.  (11)  with 


What  about  Onsager’s  regression  to  equilibrium?  We  have  to  create  an  out-of¬ 
equilibrium  condition  and  we  have  to  assess  under  which  condition  the  regression  to 
equilibrium  fits  the  prescription  of  Eq.  (4).  This  discussion  will  shed  light  into  the 
quantum-like  formalism  used  in  Section  4  to  evaluate  the  equilibrium  correlation 
function  </>*(?).  On  the  other  hand,  the  formal  approach  of  Section  4  will  make  it 
possible  for  us  to  prove  that  the  Onsager  regression  is  difficult  to  realize  in  practice. 

We  have  seen  that  the  time  evolution  of  the  correlation  function  <P^{t)  is 
determined  only  by  the  operator  T  responsible  for  the  time  evolution  of  the  anti¬ 
symmetric  part  of  the  distribution  density.  This  means  that  the  process  of  back 
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injection  taken  into  account  by  C(t),  see  Eq.  (22),  is  not  activated.  From  a  careful 
consideration  of  Eq.  (24)  the  reader  can  convince  him/herself  that  the  time  evolution 
of  pA(y,  t),  with  v<  1  is  disconnected  from  the  time  evolution  of  pA(y,  t),  with  y>  1. 
This  is  an  important  fact  that  explains  why  there  is  no  need  to  ensure  for  the  anti¬ 
symmetric  part  the  same  continuity  conditions  as  those  made  necessary  to  assign  a 
physical  meaning  to  the  symmetric  part.  In  conclusion,  we  are  allowed  to  replace 
the  formal  picture  of  Eqs.  (23)  and  (24)  with  a  reduced  picture  only  involving  the 
region  [0,1], 


Q-tPred(y>  0  =  r redPred(y,  0  , 


(42) 


where 


rred^-),l  /  .  (43) 

0y 

This  corresponds  to  the  time  evolution  of  trajectories  moving  from  a  given  initial 
position  y<  1  up  to  y  =  1,  where  the  particle  is  instantaneously  absorbed. 

To  shed  light  into  the  somewhat  disconcerting  negative  probability  distribution, 
concerning  the  right  part  of  the  anti-symmetric  probability  distribution,  we  can 
imagine  the  particles  on  the  right  as  being  holes.  The  holes  move  from  y>  1  to  the 
left,  toward  y  —  1,  while  the  particles  move  from  y<  1,  to  the  right,  toward  y  —  1, 
again.  Particles  and  holes  annihilate  one  another  at  y  =  1 .  This  particle-hole 
representation  only  refers  to  the  anti-symmetric  part  of  the  probability  distribution. 
In  this  case,  by  definition,  the  number  of  particles  is  identical  to  the  number  of  holes, 
so  making  the  conservation  of  the  number  of  particles  depend  only  on  Eq.  (22).  Eq. 
(22)  ensures  the  conservation  of  the  number  of  particles  by  means  of  the  back 
injection  process,  which  assigns  to  the  right  and  to  the  left  the  same  number  of 
particles  per  unit  of  time.  The  total  number  of  particles  and  holes  decreases  due  to 
the  mutual  annihilation  process,  thereby  yielding  the  relaxation  of  the  mean  value  of 

(£).  This  variable  gets  the  value  1  for  a  particle,  and  the  value  —1  for  a  hole.  Thus, 
(£)  is  proportional  to  the  total  number  of  particles  and  holes,  thereby  undergoing  a 
relaxation  process.  In  other  words,  rather  than  introducing  the  disturbing  concept  of 
negative  probability  distribution  for  the  anti-symmetric  part,  we  introduce  the 
notion  of  hole.  Both  particles  and  holes  are  characterized  by  positive  probability 
distributions,  with  the  particles  moving  from  the  left  to  the  right  and  the  holes  in  the 
opposite  direction.  All  this  is  formal,  and  serves  the  purpose  of  clarifying  the  twofold 
nature  of  the  equilibrium  correlation  function,  an  equilibrium  property  taking  care 
of  an  out-of-equilibrium  process  such  as  the  regression  to  equilibrium  subsequent  to 
a  perturbation. 

Let  us  first  of  all  discuss  a  physical  condition  where  the  prescription  of  Eq.  (4)  is 
fulfilled.  Let  us  consider  the  initial  distribution  p(y,  0)  defined  as  follows: 

p(y,  0)  =  2p0(y),  y  <1,  (44) 

where  pQ(y)  denotes  the  equilibrium  distribution  of  Eq.  (20) 

p(y,  0)  =  o,  y>  1  . 


(45) 
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Note  that  the  factor  of  2  is  made  necessary  by  the  fact  that  we  are  moving  the 
particles  located  on  the  right  region,  at  2  —  y,  to  the  left,  at  y,  with  y  <  1 .  This  has  the 
effect  of  creating  an  asymmetric  initial  condition,  where  the  number  of  particles  is 
not  compensated  by  the  number  of  holes.  The  anti-symmetric  part  of  this  out-of¬ 
equilibrium  initial  condition,  where  the  number  of  particles  is  identical  to  the 
number  of  holes,  coincides  with  the  anti-symmetric  part  of  the  equilibrium 
distribution.  With  this  merely  formal  trick  we  make  easier  the  problem  of  evaluating 
(£},  which  is  equal  to  the  population  difference  between  the  left  and  the  right  state, 
denoted  by  us  as  A P(t).  It  is  evident  that 


d ypiv,  t)  , 


(46) 


where  p(y,  t)  is  driven  by  the  reduced  picture  of  Eq.  (43).  It  is  also  evident  that 


exp(  r  re<i  t  )p(y,  0 )  dy  =  2 


exp(T t)pxAq\y,  0)  d y  , 


(47) 


where  p^'iy,  0)  denotes  the  left  portion  of  the  anti-symmetric  distribution 
corresponding  to  the  equilibrium  distribution.  On  the  basis  of  the  results  of  Section 
4,  we  know  that  fJ  exp(/Y)p^?)(y,  0)dy  is  the  equilibrium  correlation  function.  On 
the  other  hand,  as  pointed  out  earlier,  A P(t)  =  (£).  Of  course,  this  means  that  y<  1 
corresponds  to  ^  =  1  and  y  >  1  to  £  =  —  1 .  Thus  we  conclude  that  in  this  case  the 
Onsager  condition  of  Eq.  (4)  is  exactly  fulfilled,  with  K  —  2. 

What  about  the  regression  to  equilibrium  in  general?  For  any  initial  distribution 
p(y,  0)  we  adopt  the  same  approach  as  the  one  earlier  used.  We  create  a  new  initial 
distribution,  with  the  right  part  empty,  such  that  the  anti-symmetric  portion  of  the 
new  out-of-equilibrium  distribution  is  identical  to  the  anti-symmetric  portion  of  the 
original  distribution.  Thus,  we  prove  that 


A  Pit)  =  2  f  exp(T  t)p{nAoneq)(y,  0)  dy  .  (48) 

Jo 

Note  that  with  p^oneq'(y,  0)  we  denote  the  left  portion  of  the  anti-symmetric  part  of 
the  initial  distribution  p(y,  0)  when  this  is  not  the  equilibrium  condition.  How  to 
realize  this  condition  in  such  a  way  as  to  fulfill  the  Onsager  principle?  A  possible 
answer  can  be  found  in  Ref.  [7].  The  authors  of  this  paper  adopted  a  deterministic 
intermittent  map,  of  which  our  dynamic  model  of  Section  2  is  an  idealization.  The 
geometric  perturbation  prescription  adopted  by  these  authors  is  equivalent  to 
imagine  that  at  a  time  ta<  0  the  parameter  2  of  left  portion  is  abruptly  turned  into 
2'  ^2.  If  ta  =  — oo,  the  system  has  enough  time  to  regress  to  the  corresponding 
equilibrium.  In  this  equilibrium  condition,  the  left  and  right  distribution  have  the 
same  analytical  form,  but  different  statistical  weight.  This  can  be  described  by 


Py  =  (2  -  r) 


0(1 -y),  0(y  —  1 ) 

Pl - ; — i - b  Pr - : — r 

y  1  (2  -  y) 


(49) 
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with 

Pr+Pl  =  1  (50) 

and 

Pr^Pl  ■  (51) 

It  is  evident  that  this  special  initial  condition  fulfills  the  Onsager  principle.  In  fact  in 
this  case 

/req)(y,o)=p0(y).  (52) 

However,  for  this  property  to  hold  true  we  have  to  assume  that  the  perturbed 
condition  applies  from  t  —  ta=  —  oo  to  t  —  0.  What  about  making  ta>o o  and  closer 
and  closer  to  t  =  0?  We  argue  that  this  yields  the  breakdown  of  the  Onsager 
principle.  Note  that 

dA^(?>  =  2  J  rred  exp(T redt)Pred(y>  0)  d )’  =  pred(  1,  t )  .  (53) 

This  suggests  that  the  relaxation  time  corresponds  to  the  time  of  transit  from  y<  1  to 
y—  1,  where  the  particle  is  instantaneously  annihilated.  On  the  other  hand,  an 
averaging  process  must  be  carried  out,  which  has  to  take  into  account  the 
distribution  of  particles  at  t  —  0.  We  think  that  making  ta  finite  has  the  effect  of 
rendering  incomplete  the  evolution  toward  the  equilibrium  distribution  in  the 
presence  of  perturbation.  If  pred(y,  0)  =  p()(y),  the  regression  of  pred(  1 ,  t)  is 
proportional  to  the  equilibrium  correlation  function.  However,  if  the  regression  to 
equilibrium  in  the  presence  of  perturbation  is  incomplete,  the  time  evolution  toward 
equilibrium  must  be  evaluated  in  a  quite  different  way.  We  assume  that  in  the  regions 
close  to  the  border  y  =  0,  the  initial  distribution  deviates  from  the  equilibrium 
prescription  of  being  proportional  to  l/yz~l.  More  precisely,  we  assume  that  for 
e<y<  1,  the  distribution  p(y,t )  is  proportional  to  l/yz~x,  whereas  for  the 
distribution  might  or  might  not  depart  from  equilibrium.  In  the  case  where  it  does, 
the  distribution  can  be,  for  instance,  independent  of  y.  In  the  ideal  case  ta  —  — oo, 
s  —  0,  with  no  breakdown  of  the  Onsager  principle. 

To  study  the  breakdown  of  this  principle,  we  evaluate  the  asymptotic  behavior  of 
A P(t).  We  select  a  set  of  initial  conditions  j(0)’s  with  O^jfO)^;:.  We  call  M  the 
number  of  particles  located  in  this  region  at  t  =  0.  The  asymptotic  behavior  of  A  P(t) 
is  determined  by  the  time  necessary  for  these  trajectories  to  reach  the  border  y  =  1 . 
The  first  trajectory  will  reach  the  border  after  a  given  time  t  —  T first,  after  which  the 
number  M  will  begin  decreasing,  thereby  determining  the  decay  of  A P(  t)  in  this  time- 
asymptotic  region.  It  is  straightforward  to  prove  that  the  time  of  arrival  at  y  —  1  of 
the  trajectory  with  initial  condition  y(0),  called  t,  is  related  to  j(0)  by 


y( 0; t ) 


1 

[1  +(z-  ' 


(54) 


Thus  we  obtain 


d M  dy(0;  t) 
d  t  d? 


t0’(0;  0)  • 


(55) 
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First  of  all,  let  us  show  that  this  procedure  recovers  the  result  of  the  earlier  two 
sections.  Let  us  assume  that  p(y(0;t))  is  always  proportional  to  I  />'(0;  f)z_1 .  In  this 
case  we  obtain 

AM  1 

—  = - r  .  (56) 

At  [1  +  (z  -  1 

Then  we  evaluate  M(t)  from 

r'  AM 

M{t)=l+l  d?V  •  (57) 

By  using  Eq.  (57)  we  recover  the  Onsager  principle,  with  M(t)  oc  to  l/d1-2.  In  the 
case  where  p(y( 0;  t))  is  independent  of  j(0;  ?),  the  same  calculation  procedure  yields 
M(t)  oc  1  If--1,  which  sanctions  the  breakdown  of  the  Onsager  principle. 


6.  Higher-order  correlation  functions 


We  now  show  that  the  calculation  of  higher  order  correlation  functions,  though 
difficult,  can  be  done  by  establishing  an  even  deeper  connection  with  the  quantum 
mechanical  perspective.  Let  us  address  the  problem  of  evaluating  the  fourth-order 
correlation  function  ■  As  shown  in  Section  4,  we  move  from 

the  equilibrium  distribution  and  let  it  evolve  for  a  time  t\ .  The  distribution  selected 
is  in  equilibrium.  Therefore  it  will  remain  unchanged.  At  time  t  —  t\  we  apply 
the  operator  £  to  the  distribution.  Since  the  equilibrium  distribution  is  symmetric, 
the  application  of  the  sign  operator  changes  it  into  the  anti-symmetric  dis¬ 
tribution.  We  let  this  distribution  evolve  for  the  time  tj  —  t\.  This  has  the  effect  of 
yielding 

A)+ti°neq)(y,t2-h). 

At  this  stage  there  are  two  possibilities: 


(a) 

(b)  pfne(l) 


(y,t2  -  h)  =  0; 
(y,t2  -  1 1)#0. 


Let  us  consider  case  (a)  first.  In  this  case,  we  proceed  as  follows.  At  time  L  we 
apply  to  the  distribution  the  operator  £,  and  we  change  it  into  the  original 
equilibrium  distribution.  This  means  that  the  time  evolution  from  L  to  h  leaves  it 
unchanged.  At  time  we  apply  to  it  the  operator  C,  and  we  turn  it  into  p^q)  again. 
We  let  this  distribution  evolve  till  time  t*.  At  this  time  we  apply  to  it  the  operator 
again  and  make  the  final  average.  The  result  of  condition  (a)  yields 


(mmmm)  =  m)m)mw  i)> . 


(58) 
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However,  condition  (a)  applies  only  when  z  =  1.  In  this  case  it  is  straightforward  to 
prove  that  the  equilibrium  distribution  is  flat,  the  anti-symmetric  component  consists 
of  two  flat  components,  positive,  y>  1,  and  negative,  y<  1,  which  are  left  unchanged 
by  the  application  of  the  operator  T  of  Eq.  (25). 

In  the  case  r>  1  condition  (a)  is  violated  and  with  it  the  factorization  condition  of 
Eq.  (58).  This  agrees  with  the  recent  work  (21),  where  an  explicit  expression  for  the 
fourth-order  correlation  function  in  terms  of  $=(?)  was  given  by  means  of  an 
approximated  approach  that  nevertheless  leads  to  an  accurate  reproduction  of 
numerical  results.  The  demonstration  was  made  by  applying  the  method  of 
conditional  probabilities  to  the  study  of  single  trajectories.  Thus  we  are  forced  to 
consider  case  (b). 

The  problem  with  condition  (b)  is  that  it  yields  a  distribution  with  the  symmetric 
component  departing  from  equilibrium.  This  departure  is  in  an  apparent  conflict 
with  the  assumption  that  the  correlation  function  is  calculated  using  the  equilibrium 
condition.  In  fact  it  seems  to  be  equivalent  to  stating  that  the  calculation  of  an 
equilibrium  correlation  function  generates  an  out-of-equilibrium  condition.  Let  us 
see  why  this  appears  to  be  true. 

In  the  case  of  the  two-time  correlation  function  we  apply  the  operator  £  to  the 
density  distribution  twice.  The  first  application  allows  us  to  observe  the  time 
evolution  of  the  anti- symmetric  component  distribution  density,  with  no  conflict 
with  equilibrium,  given  the  fact  that  the  symmetric  component  remains  at 
equilibrium.  The  second  application  of  <L  turns 

-  ti)p^q\y,  0)  +p(™"eq)(y,  h  -  ti) 

into 

(l’c(t2  ~  ti)p(gq\y,  0)  +Ps°neq)(y,  t2  -  t\ )  . 

The  calculation  done  in  Section  4  proves  that  T \'\j>^"neq)(y,  t2  —  L  )]  =  0,  with  the 
symbol  Tr  denoting,  for  simplicity,  the  integration  over  v  from  0  to  2.  To  evaluate 
the  fourth-order  correlation  function,  after  applying  the  operator  £  for  the  second 
time,  we  must  study  the  time  evolution  of  p™neq\y,  t2  —  t\)  from  t2  to  t2,  yielding 
p\ ™neq )(y,  f3,  t2,  t2  —  /i).  This  is  a  contribution  generating  some  concern,  since  it 
activates  again  the  back  injection  process,  which  we  have  seen  to  be  intimately 
related  to  the  deviation  from  equilibrium.  However,  the  compatibility  with 
equilibrium  condition  is  ensured  by  the  property  Tr \p^oneq\y,  t2,  t2,  t2  —  t\ )]  =  0. 

At  time  t2  we  have  to  apply  the  operator  £  again,  and  this  allows  us  to  make  an 
excursion  in  the  anti-symmetric  representation,  with  the  time  evolution  given  by  the 
operator  expiry  —  ?3)).  At  time  ?4  we  apply  the  operator  £  again,  we  go  back  to  the 
symmetric  representation  and  we  conclude  the  calculation  by  means  of  trace 
operation. 

In  conclusion,  the  compatibility  with  equilibrium  is  guaranteed  by  the  fact  that  at 
the  intermediate  steps  of  the  calculation  Tr[;;]  =  0  (the  final  step,  of  course,  generates 
the  density-generated  correlation  function,  thereby  implying  Tr[;;]^0).  If  the 
intermediate  p  is  anti-symmetric,  this  condition  is  obvious.  If  the  intermediate  p  is 
symmetric,  the  vanishing  trace  condition  generates  the  apparently  unphysical 
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property  that  the  symmetric  contribution  gets  negative  values  over  some  portions  of 
the  interval  I.  We  have  to  stress,  however,  that  the  Liouville-like  approach  illustrated 
in  this  paper  keeps  the  distribution  density  p(y,  t)  positive  definite,  as  it  must.  The 
generation  of  a  negative  distribution  density  refers  to  the  calculation  of  equilibrium 
correlation  functions,  of  any  order,  and  it  is  a  quantum-like  property  that  must  be 
adopted  to  guarantee  that  the  genuine  distribution  density  remains  in  the 
equilibrium  condition. 

It  is  worth  pointing  out  that,  although  these  arguments  show  that  it  is  possible  to 
evaluate  the  fourth-order  correlation  function,  we  could  not  find  its  explicit 
expression,  and  as  a  consequence  it  is  not  yet  possible  for  us  to  assess  whether  or  not 
the  explicit  result  of  Ref.  [21]  can  be  recovered.  Furthermore,  the  explicit  result  of 
Ref.  [21]  is  not  exact.  Consequently,  we  have  to  conclude  that  the  evaluation  of  the 
fourth-order  correlation  function  in  the  non-Poisson  case  is  still  an  open  problem. 


7.  Comments  on  the  quantum-like  formalism 

The  reader  might  be  disconcerted  by  the  properties  of  the  quantum-like  formalism 
adopted  in  this  paper,  and  especially  by  those  illustrated  in  Section  6,  such  as  Tr[/>]  = 
0,  implying  that  also  the  symmetric  part  of  the  probability  distribution  might  become 
negative,  as  well  as  the  anti-symmetric  part,  as  discussed  in  Section  4.  Actually,  we 
did  not  derive  these  properties  from  the  Liouville-like  discussed  in  this  paper.  We 
have  been  rather  forced  to  make  the  plausible  conjecture  that  these  properties  are 
true  by  the  wish  of  extending  the  trajectory-density  equivalence  from  Poisson  to  non- 
Poisson  statistics. 

The  authors  of  Ref.  [16]  questioned  the  equivalence  between  density  and  trajectory 
picture  on  the  basis  of  a  generalized  diffusion  equation  (GDE)  derived  from  Eq.  (14), 
thought  of  as  being  the  generator  of  diffusion.  The  GDE  of  Ref.  [16]  yields  a 
dramatic  conflict  with  the  adoption  of  the  trajectory  perspective.  The  reason  for  this 
conflict  is  that  the  authors  of  Ref.  [16]  adopted  condition  (a)  of  Section  6  implicitly. 
If  we  adopt  condition  (b)  of  Section  6  instead,  to  avoid  the  trajectory-density  conflict 
emerging  from  (a),  then  we  have  to  accept  all  the  disconcerting  properties  of  the 
distribution  density  formalism.  In  Section  4  we  have  adopted  the  particle-hole 
perspective  to  justify  the  anti-symmetric  part  of  the  distribution  density  becoming 
negative.  It  is  not  enough.  The  correct  evaluation  of  higher-order  correlation 
functions  implies  that  also  the  symmetric  component  must  be  non-positive  definite, 
although  in  the  virtual  sense  discussed  in  Section  6. 

We  think  that  all  this  might  help  to  understand,  by  analogy,  the  strange  and 
apparently  paradoxical  properties  of  quantum  mechanics,  whose  statistical  treat¬ 
ments  is  based  on  the  quantum  Liouville  equation.  The  counterpart  of  this  quantum 
treatment  is  given  here  by  Eq.  (14).  To  establish  a  satisfactory  connection  between 
Liouville  and  CTRW  formalism,  we  have  to  prove  the  plausible  conjectures  of 
Section  6,  in  spite  of  their  disconcerting  nature.  The  only  possible  alternative  would 
be  to  abandon  the  Liouville  approach  (and  quantum  mechanics  with  it?)  and  to 
adopt  the  CTRW  formalism  instead. 


286 


P.  Allegrini  et  al.  /  Physica  A  347  (2005)  268-288 


8.  Concluding  remarks 

In  conclusion,  this  paper  shows  how  to  derive  the  equilibrium  correlation  function 
using  only  information  afforded  by  the  Liouville-like  approach.  The  intriguing 
problem  to  solve  was  how  to  use  the  Liouville  equation,  without  conflicting  with  the 
equilibrium  condition.  The  solution  to  this  intriguing  problem  is  obtained  by 
splitting  the  Liouville  equation  into  two  independent  components,  the  symmetric 
component  corresponding  to  Eq.  (22)  and  the  anti-symmetric  component 
corresponding  to  Eq.  (23).  This  splitting  allows  us  to  study  the  regression  to 
equilibrium  of  the  correlation  function  (Pe(t),  through  Eq.  (23),  without  ever 
departing  from  equilibrium,  a  physical  condition  that  is  controlled  by  the 
independent  equation  of  motion  for  the  symmetric  component,  Eq.  (22).  This  is 
very  formal,  and  Section  5  makes  it  possible  for  us  to  reconcile  it  with  physical 
intuition.  We  imagine  that  the  Liouville  equation  is  used  to  evaluate  the  difference 
between  the  population  of  the  left  and  right  state.  This  makes  it  possible  to  establish 
a  direct  connection  between  the  experiment  of  regression  to  equilibrium  and  the 
formalism  of  Sections  2  and  4.  We  have  to  create  an  asymmetrical  initial  condition, 
with  more  population  on  the  left  than  on  the  right.  The  time  evolution  of  A P(t) 
depends  only  on  the  time  evolution  of  the  anti-symmetric  component  of  the 
distribution  density,  and  consequently  only  on  the  operator  T.  This  is  the  reason  why 
it  is  possible,  in  principle,  to  connect  regression  from  an  out-of-equilibrium  initial 
condition  to  the  equilibrium  correlation  function,  a  condition  that  implies  no 
deviation  from  equilibrium  distribution.  Elowever,  creating  in  a  finite  time  an  out-of¬ 
equilibrium  condition  such  that  the  left  part  of  the  anti-symmetric  component  is 
identical  to  the  left  part  of  the  equilibrium  distribution  is  impossible.  This  is  the 
reason  why  we  predict  the  breakdown  of  the  Onsager  principle  in  general. 

Section  6  explains  why  in  the  literature  on  dichotomous  fluctuations  the 
factorization  assumption  of  Eq.  (58)  is  often  made  regardless  of  the  Poisson  or 
non-Poisson  nature  of  the  underlying  process.  See,  for  instance,  the  work  of  Fulinski 
[25]  as  well  as  Ref.  [16].  In  fact,  if  condition  (a)  applies,  the  higher-order  correlation 
functions  are  factorized,  thereby  making  their  calculation  easy.  However,  this 
assumption  conflicts  with  the  trajectory  arguments  of  Ref.  [21],  which  prove  that  the 
factorization  condition  is  violated  by  the  non-Poisson  condition.  A  rigorous  use  of 
the  Liouville  equation  shows  that  condition  (a)  does  not  apply,  and  that  we  have  to 
use  condition  (b)  instead.  The  calculation  of  the  fourth-order  correlation  function  is 
not  straightforward,  and  this  is  the  reason  why,  to  the  best  of  our  knowledge,  it  was 
never  done  using  the  Liouville  approach.  Furthermore,  as  earlier  pointed  out  in 
Section  7,  the  rules  discussed  in  Section  6  are  conjectured  rather  than  proved,  and 
correspond  to  the  mathematical  properties  to  guarantee  the  density-trajectory 
equivalence. 

This  is  a  fact  of  some  relevance  for  the  creation  of  master  equations  with  memory. 
There  are  two  major  classes  of  GME.  The  first  class  is  discussed,  for  instance,  in  Ref. 
[1].  The  master  equations  of  this  class  are  equivalent  to  the  CTRW  of  Montroll  and 
Weiss  [3]  and  are  based  on  the  waiting  time  distribution  The  second  class  of 
master  equations  is  based  on  the  correlation  function  ^(f).  Recent  examples  of  this 
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second  class  can  be  found  in  Ref.  [16]  and  in  Ref.  [26].  Due  to  the  direct  dependence 
on  the  correlation  function  <IJc(t),  the  derivation  of  the  master  equations  of  this 
second  class  is  made  easy  by  the  factorization  assumption.  It  must  be  pointed  out,  on 
the  other  hand,  that  the  factorization  property,  which  is  not  legitimate  with  renewal 
non-Poisson  processes,  is  a  correct  property  if  the  deviation  from  the  exponential 
relaxation  is  obtained  by  time  modulation  of  a  Poisson  process  [26].  Beck  [27]  is  the 
advocate  of  the  modulation  process  as  generator  of  complexity.  Thus,  we  find  that 
the  master  equations  of  the  first  class  are  generated  by  the  renewal  perspective  of 
Montroll  and  Weiss  [3]  and  those  of  the  second  class  are  the  appropriate  tool  to 
study  complexity  along  the  lines  advocated  by  Beck  [27].  We  think  that  the  results  of 
the  present  paper  might  help  the  investigators  in  the  field  of  complexity  to  make  the 
proper  choice,  either  modulation  or  renewal  [28],  or  a  mixture  of  the  two  conditions. 

As  a  final  remark,  let  us  go  back  to  the  issue  raised  by  Ref.  [16]  about  the 
trajectory-density  conflict.  On  the  basis  of  the  results  of  this  paper  we  can  conclude 
that,  at  first  sight,  this  conflict  can  be  judged  to  be  apparent,  insofar  as  the 
conclusions  of  the  authors  of  this  paper  were  based  on  the  assumption  that  the 
Liouville-like  approach  yields  the  factorization  property  of  Eq.  (58),  which  is  correct 
only  in  the  Poisson  case.  However,  it  is  not  yet  possible  to  express  a  final  verdict  on 
this  important  issue.  This  is  so  because  the  trajectory  method  [21]  does  not  yield  yet 
an  exact  expression  for  the  fourth-order  correlation  function,  and  using  the 
arguments  of  Section  6  we  could  not  arrive  yet  at  any  expression,  either  exact  or 
approximated,  from  the  density  perspective.  In  Section  6  we  have  seen  that  the 
calculation  of  the  fourth-order  correlation  function  activates  the  contribution  C(t), 
which  is  responsible  for  an  aging  effect  [1].  On  the  other  hand  the  surprising  effects 
[28]  discovered  by  Sokolov,  Blumen  and  Klafter  [12]  indicate  that  the  non-Poisson 
property-induced  aging  effect  might  lead  us  to  adopt  with  caution  the  Liouville-like 
approach.  Similar  conclusions  have  been  recently  reached  by  Aquino  and  coworkers 
[30]  who  have  been  forced  to  study  the  absorption  process  of  radiation  by  a  non- 
Poisson  system  with  a  method  entirely  based  on  the  CTRW,  given  the  fact  that  the 
ordinary  Liouville-like  approach  cannot  properly  take  into  account  the  aging  effect. 
In  conclusion,  we  cannot  rule  out  the  possibility  that  the  issue  of  the  trajectory- 
density  conflict  raised  by  Ref.  [16]  might  signal  the  symptoms  of  a  disease  that  can 
only  be  cured  by  moving  from  the  Liouville-like  to  the  CTRW  perspective. 


Acknowledgements 

PG  thankfully  acknowledges  the  Army  Research  Office  for  the  financial  support 
through  Grant  DAAD 19-02-0037. 


References 

[1]  P.  Allegrini,  G.  Aquino,  P.  Grigolini,  L.  Palatella,  A.  Rosa,  Phys.  Rev.  E  68  (2003)  056123. 

[2]  R.  Kubo,  M.  Toda,  N.  Hashitsume,  Statistical  Physics  II.  Nonequilibrium  Statistical  Mechanics, 
Springer,  Berlin,  1991. 


288 


P.  Allegrini  et  al.  /  Pliysica  A  347  (2005)  268-288 


[3]  E.W.  Montroll,  G.H.  Weiss,  J.  Math.  Phys.  6  (1965)  167. 

[4]  J.  Klafter,  M.  Shlesinger,  G.  Zumofen,  Phys.  Today  49  (2)  (1996)  33. 

[5]  R.  Metzler,  J.  Klafter.  Phys.  Rep.  339  (2000)  1. 

[6]  B.B.  Mandelbrot,  The  Fractal  Geometry  of  Nature,  W.H.  Freeman  and  Company,  New  York,  1983. 

[7]  G.  Trefan,  E.  Floriani,  B.J.  West.  P.  Grigolini.  Phys.  Rev.  E  50  (1994)  2564. 

[8]  E.  Barkai,  J.  Klafter,  Phys.  Rev.  Lett.  79  (1997)  2245. 

[9]  E.  Barkai,  J.  Klafter,  Phys.  Rev.  E  57  (1998)  5237. 

[10]  E.  Barkai,  V.N.  Fleurov,  Phys.  Rev.  E  58  (1998)  1296. 

[11]  R.  Metzler,  J.  Klafter,  I.M.  Sokolov,  Phys.  Rev.  E  58  (1998)  1621. 

[12]  I.M.  Sokolov,  A.  Blumen,  J.  Klafter,  Europhys.  Lett.  56  (2001)  175. 

[13]  J.P.  Bouchaud.  A.  Georges,  Phys.  Rep.  195  (1990)  127. 

[14]  E.  Barkai,  J.  Klafter,  Phys.  Rev.  Lett.  81  (1998)  1134. 

[15]  F.  Amblard,  A.C.  Maggs,  B.  Yurke,  A.N.  Pargellis,  S.  Leibler,  Phys.  Rev.  Lett.  81  (1998)  1135. 

[16]  M.  Bologna,  P.  Grigolini,  B.J.  West,  Chem.  Phys.  284  (2002)  115-128. 

[17]  G.  Zumofen,  J.  Klafter,  Phys.  Rev.  E  47  (1993)  851. 

[18]  T.  Geisel,  J.  Nierwetberg,  A.  Zacherl,  Phys.  Rev.  Lett.  54  (1985)  616. 

[19]  D.J.  Driebe,  Fully  Chaotic  Maps  and  Broken  Time  Symmetry,  Kluwer  Academic  Publishers, 
Dordrecht,  1999. 

[20]  X.  Brokman,  J.-P.  Hermier,  G.  Messin,  P.  Desbiolles,  J.-P.  Bouchaud,  M.  Dahan,  Phys.  Rev.  Lett.  90 
(2003)  120601. 

[21]  P.  Allegrini,  P.  Grigolini,  L.  Palatella,  B.J.  West,  Phys.  Rev.  E.,  accepted  for  publication. 

[22]  M.  Bologna,  P.  Grigolini,  M.  Karagiorgis,  A.  Rosa,  Phys.  Rev.  E  64  (2001)  016223. 

[23]  M.  Ignaccolo.  P.  Grigolini,  A.  Rosa,  Phys.  Rev.  E  64  (2001)  026210. 

[24]  N.  Goldenfeld,  Lectures  on  Phase  Transitions  and  the  Renormalization  Group,  Addison-Wesley 
Publisher,  Cambridge,  1992. 

[25]  A.  Fulinski,  Phys.  Rev.  E  52  (1995)  4523. 

[26]  M.  Bologna,  P.  Grigolini,  M.G.  Pala.  L.  Palatella,  Chaos,  Solitons  and  Fractals  17  (2003)  601. 

[27]  C.  Beck,  Phys.  Rev.  Lett.  87  (2001)  180601. 

[28]  In  the  case  where  the  renewal  condition  applies,  we  want  to  point  out  that,  as  shown  in  (1),  the  GME 
can  be  derived  from  the  CTRW  [3],  It  is  still  unknown  how  to  do  that  using  the  Liouville-like 
approach  of  Section  2,  probably  due  to  the  complexity  of  the  problem  illustrated  in  Section  6. 
Furthermore,  the  authors  of  Ref.  [12]  made  the  following  impressive  discovery.  In  the  non-Poisson 
case  even  if  a  correct  GME  is  available  [29],  its  response  to  external  perturbations  departs  from  the 
correct  prediction  that  one  can  obtain  by  perturbing  the  CTRW  instead.  Thus  the  trajectory-density 
conflict  [16]  seems  to  be  the  manifestation  of  technical  problems,  emerging  from  the  density 
perspective,  which  can  be  settled  using  the  trajectory  (CTRW)  perspective. 

[29]  Within  the  theoretical  framework  of  this  paper,  the  fractional  Fokker-Planck  equation  used  by 
Sokolov  et  al.  [12]  belongs  to  the  class  of  GDE  that  can  be  derived  from  the  CTRW. 

[30]  G.  Aquino,  L.  Palatella.  P.  Grigolini.  Phys.  Rev.  Lett.  93  (2004)  50601. 


